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Abstract. Highly accurate measurements of quantum level energies in molecular systems 
provide a test ground for new physics, as such effects could manifest themselves as minute shifts 
in the quantum level structures of atoms and molecules. For the lightest molecular systems, 
neutral molecular hydrogen (Hi, HD and D 2 ) and the molecular hydrogen ions (Hif, HD + and 
D,|), weak force effects are several orders weaker than current experimental and theoretical 
results, while contributions of Newtonian gravity and the strong force at the characteristic 
molecular distance scale of 1 A can be safely neglected. Comparisons between experiment and 
QED calculations for these molecular systems can be interpreted in terms of probing large 
extra space dimensions, under which gravity could become much stronger than in ordinary 
3-D space. Under this assumption, using the spectra of Hi we have derived constraints on 
the compactification scales for extra dimensions within the Arkani-Hamed-Dimopoulos-Dvali 
(ADD) framework, and constraints on the brane separation and bulk curvature within the 
Randall-Sundrum (RS-I and RS-II) frameworks. 


1. Introduction 

This presentation is based on work done in collaboration with Bert Schellekens, from the Institute 
of Subatomic Physics Nikhef (Amsterdam), together with Edcel Salumbides and Wim Ubachs, 
from the Vrije Universiteit of Amsterdam. In the present study, accurate results from precision 
molecular spectroscopy are exploited to put constraints on some higher-dimensional theories. 
The idea is that effects of new physics could manifest themselves as minute shifts in the quantum 
level structures of atoms and molecules. To detect such phenomena, spectroscopic experiments 
must be performed at extreme precision, where key ingredients are ultra-stable lasers and laser 
metrology techniques, and comparisons must be made on systems for which very accurate first- 
principles calculations can be performed. 

In recent years great progress has been made on the calculation of level energies in the 
neutral hydrogenic molecules (H2, HD and D2) and their ions (H((, HD + and D^"). The 
results from the variety of experimental precision measurements on both the ionic and neutral 
hydrogenic molecules are in perfect agreement with the QED-calculations, within combined 
uncertainty limits from theory and experiment. This agreement between experiment and first- 
principles calculations on the level energies of hydrogenic molecules can be used to constrain 
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the parameters associated with Physics Beyond the Standard Model, such as possible fifth 
forces between nucleons or the existence of extra space dimensions. Since weak, strong and 
(Newtonian) gravitational forces have negligible contributions to the molecular quantum level 
structure, electromagnetism is the sole force acting between the charged particles within light 
molecules. This makes it possible to derive bounds on possible fifth forces and extra dimensions 
from molecular precision experiments compared with QED-calculations El El. 

Although in the framework of String Theory extra space dimensions were found to be 
necessary from the very begining in the middle seventies of last century, it took much longer to 
appreciate the beauty and usefulness of extra dimensions in ordinary Quantum Field Theory. 
Only in the late nineties, after the seminal works of Antoniadis, Arkani-Hamed, Dimopoulos, 
and Dvali (ADD model) [3] IS] [6] on one side, and Randall and Sundrum (RS-I and RS-II 
models) [71 [8j on the other, the particle physics community started considering extra dimensions 
as a real scientific possibility deserving further investigation. The mathematical formalisms of 
these models can be applied to molecular physics test bodies, from which constraints on the 
compactification scales can be derived for the ADD scenario, whereas constraints on the brane 
separation and bulk curvature can be derived for the RS-I and RS-II scenarios. 


2. The ADD Scenario 

In the model put forward in 1998 by Arkani-Hamed, Dimopoulos, and Dvali [3, [5>J (plus 
Antoniadis [6] in its implementation in string theory) our observable universe consists of a 
3-brane embedded in a higher-dimensional cosmos. There are n extra space dimensions where 
only gravity can escape as the particles and interactions of the Standard Model remain confined 
in the ordinary 3-D space (up to a small, negligible penetration inside the bulk). The extra 
dimensions are flat and compactified with the same radius R , just for simplicity. Therefore the 
extradimensional space is an n-dimensional torus where gravitation propagates rendering the 
observed gravitational strength much weaker as the flux lines dilute into the higher-dimensional 
volume. Moreover, the propagation of gravity across the extradimensional space results in a 
modification of Newton’s law. The exact ADD gravitational potential reads [5]: 

R ADD (r) = - GY, e (2 " |fc|/R)r 

( ki,k2,...k „) 

It has the contribution from the usual massless gravitons plus a tower of Yukawa potentials 
mediated by all the massive modes: the Kaluza-Klein higher-dimensional gravitons, with 
momenta in the extra dimensions quantized in units of 2i t/R, and therefore with mass splittings 
of 27 t/R from our 3-D space perspective (from the higher-dimensional perspective they are 
massless). For distances much larger than the compactification length range, r R, the 
Newtonian potential is recovered: 
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where Mp\ is the Planck mass (10 16 TeV). For distances of the same size or shorter than the 
compactification length range, r < R, the ADD potential takes the simple form: 
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where M ( 4An ) is the fundamental mass in the (4+7i)-dimensional spacetime. Observe that this 
expression for the gravitational potential can be interpreted in the context of the Gauss law, 
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Figure 1. (Color online) Wave functions for H 2 in the electronic ground state with J = 0 and 
vibrational levels v = 0,1. 


without the need to invoke quantization of gravity. The Planck mass Mp\ is then related to the 
fundamental higher-dimensional mass Ai^ +n ^ via: 

Mk = M^ n) {R) n . (4) 

Thus the fundamental mass M ( 4+n ) may still be small and Mp\ becomes large due to the 
compactified volume of extra dimensions. Arkani-Hamed et al. proposed that Mr 4 +n ) equals 
the electroweak scale Mew hi order to solve the hierarchy problem, explaining in this manner 
the desert between the electroweak scale of 1 TeV and the Planck scale of 10 16 TeV. In the 
present study, since it is not our intention to solve the hierarchy problem, we leave Mu +n \ as 
a free parameter, together with R and n, and we attempt to constrain these parameters from 
molecular physics experiments. 

The gravitational potential has an effect on the level energy of a molecular quantum state 
with wave function \E r (r), written as an expectation value [lj: 
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'P(r) represents the probability that the internuclear separation is r, the vibrational 
coordinate that probes the forces between the nucleons. The first integral represents the ordinary 
gravitational attraction, which for protons is 8 X 10~ 37 times weaker than the electrostatic 
repulsion, and can therefore be neglected. The second integral represents the effect of modified 
gravity and is evaluated using accurate wave functions for the molecules. The wave functions 
of the H 2 ground electronic state for J = 0 and vibrational levels v = 0,1 are plotted in Fig. [l] 
as obtained from ab initio calculations mm ■ In practice, the integration is performed up to 
r = 10 A since the wave function amplitude is negligible beyond that. In addition, at distances 
r < 0.1 A the wave function amplitude becomes negligible so that the second integral in Eq. ([5J) 
converges without additional assumptions. 

For spectroscopic transitions between quantum states Ti and T 2 in molecules, the expectation 
value for a higher-dimensional gravity contribution is given by a differential effect: 
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Here the ADD-expectation value is written explicitly as a function of the two relevant 
parameters: the number n of extra space dimensions and the compactification scale R. 
Obviously, measurements on the dissociation limit of molecules, where 4/ i is the lowest energy 
bound state and 4/2 is the non-interacting two-atom limit at r = 00 , are the most sensitive 
probes as the diference between 4 r and 4/2 is larger. 


3. The Randall-Sundrum Scenarios 

Let us now consider the Randall-Sundrum scenarios, RS-I [Tj and RS-II [8], presented by the 
authors in 1999. They consist of two parallel 3-branes separated by a distance y c along one extra 
dimension y. The branes have opposite tension T (vacuum energy) and the bulk has a negative 
cosmological constant. According to General Relativity, where all types of energy are gravity 
sources, the branes and the bulk produce gravitational fields, giving rise to the AdS metric: 

ds 2 = e- 2k \ y \r hlv dx tl dx v + dy 2 , (7) 


where exp(— k\y\) is the warp factor and k is the bulk curvature. The strength of gravity is 
highest on the positive T > 0 tension brane and falls off along the extra dimension. 

In analogy with the ADD model, in the RS models the particles and interactions of the 
Standard Model are confined in one 3-brane, that we denote the SM-brane, whereas gravity 
propagates everywhere. Observe that the warped metric crucially differentiates the RS models 
from the ADD model with a flat metric where k = 0. 

The modified gravitational potential between two masses separated a distance r in the SM- 
brane can be expressed as: 

^Rs(r) = -G rn ^ 2 (l + A RS ), (8) 

where Ars is the RS-I or RS-II correction to the Newtonian potential. 

In the RS-I scenario the SM-brane is the negative T < 0 tension brane. The authors did not 
present any gravitational potential for this model, however, and it was Callin [9], in 2004, who 
presented its computation. For short distances, ir<l, he obtained: 
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where one can distinguish two regimes, ky c <C 1 and ky c 3> 1, with the result at leading orders: 
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The potential for long distances, lr> 1, does not apply for the scales probed by molecules (it 
gives inconsistent results) and we do not consider it further. 

In the RS-II scenario the SM-brane is the positive T > 0 tension brane. The T < 0 tension 
brane can be set to infinity, y c —> 00 , as the massless gravitons are localized very close to the 
T > 0 brane, their wave function decreasing drastically along the extra dimension y. Therefore, 
in the RS-II scenario the extra dimension needs not be compactified. The authors derived a 
potential with the long distance limit that has been applied to astrophysics (see for example mi). 
For short distances, kr <C 1, Callin and Ravndal [TO] obtained: 
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Note the correspondence of this result with that of Eq. (10) in the limit ky c 1. This limit 
together with h<l seems to indicate a smooth transition from RS-I to RS-II for small bulk 
curvature k <C 1/r and y c 3> l/k, as the brane separation y c is increased to infinity [TO]. 

From the RS corrections to the Newtonian potential the expectation values of the leading- 
order shifts of transitions in molecules, for short distances kr <C 1, are therefore: 
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where T = (1 — e~ 2kyc )/(l + | e ~ 2kyc ) for the RS-I scenario and T = 1 for the RS-II scenario. 
Using these expressions, we derive bounds on the brane separation y c and the bulk curvature k 
from molecular spectroscopy data. 

The RS-II scenario with a single brane can be generalized to n extra dimensions provided 
that the mixing effects of the warp factors for each extra dimension are small compared with 
the leading order terms. In this case the resulting correction to the Newtonian potential is, to 
a good approximation, the product of the leading corrections of each extra dimension. For n 
extra dimensions, if the bulk curvatures have the same value k, the RS-II correction is then 
Ars-ii,™ = (ARg_n) n , where we use Ars-ii in Eq. 0- 


4. Constraints on extra dimensions from molecular spectroscopy data 

In the previous sections, the expectation value for a higher-dimensional gravity contribution 
to a transition frequency in a molecule was presented for both ADD and RS approaches to 
higher dimensions. This expectation value is interpreted as a contribution to the binding 
energy of molecules in certain quantum states. However, as discussed in Ref. [Tj , analysis of 
recent experiments with hydrogen neutral molecules and their ions, together with their stable 
isotopomers containing deuterons, showed that QED computations are in very good agreement 
with observations for these molecular systems. Consequently, one can derive constraints on the 
parameters underlying the higher-dimensional theories: the compactification scale R depending 
on the number of extra dimensions n, for the ADD scenario, and the bulk curvature k or brane 
separation y c for the RS scenarios. For that purpose one can use the combined uncertainty 5E 
for molecular systems, where <5F t h e ory and SE exp are uncertainties of theory and experiment: 

^ = sJSEl xp + 5E 2 heory . (13) 

In the ADD scenario the constraining relation 
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translates on constraints on the compactification scale R, depending on n, via: 
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where A is the difference in expectation values over the wave function densities in the molecule: 



(16) 


As shown in Fig. [lj the wave functions for the lowest vibrational states, in the case of 
H 2 and for J = 0, are located in the same region of space. Consequently, the fundamental 
vibrational transition in the hydrogen molecule (v = 0 —> v = 1) probes only a differential effect. 
Constraints on R, obtained from a comparison with the fundamental vibrational transition of 
H 2 are presented in Table [l] 




















Table 1. Constraints on the compactification scale R of extra dimensions (in units of m) as 
derived from both the measurement of the fundamental (v = 0 —> 1) vibration in H2 and of 
the dissociation limit D 0 of H2. The constraints are derived within the ADD-scenario assuming 
that the n extra dimensions are of equal size. The higher-dimensional Planck length Rpiy4+ n ) 
(in units of m) and Planck mass M( 4+ri ) (in units of GeV) for Dq are also tabulated. 
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The experimental as well as the theoretical results for the fundamental vibration in the 
hydrogen molecule are known to the 10” 4 cnr 1 level, an order of magnitude more accurate than 
the values for the binding energies (181 [19i l. However for a comparison of dissociation limits, 
it is no longer a small difference along the internuclear coordinate axis that is probed, but the 
difference between the molecular scale of 1 A and infinite atomic separation. The expectation 
value for the ADD-contribution to the binding energy of the lowest bound state in the H2 
molecule, the Dq binding energy, is: 

(AP A DD(n, R)) = Gmim 2 R n (17) 

\ r / fo 

By comparing to the experimental findings on Ho(H2) [16j this leads to another set of constraints 
on R for n extra dimensions, which are also listed in Table [I] In Ref. [T] more constraints on 
R are presented corresponding to the dissociation limits of the deuterium molecule D 2 and the 
(v = 0 —> 4) R(2) ro-vibrational transition in HD + . 

Similarly, one can derive constraints pertaining to corrections in the RS models and the 
combined uncertainty 5E for a specific molecular transition 


(A V RS (k)) < 5E. 


(18) 


Using 5E for the dissociation limit Ho(H2) we obtain the following constraints. For the RS-I 


scenario in the short distance (kr <C 1) regime, using Eq. (10) one finds a constraint on the 
brane separation of y c < 2 X 10 18 m in the limit ky c <C 1, and a constraint for the inverse of the 
curvature of 1/k < 3 x 10 18 m in the limit ky c 1. 

For the RS-II scenario with one extra dimension in the short distance (kr <C 1) regime, using 
Eq. (12) one finds a constraint for the inverse of the curvature of 1/k < 3 X 10 18 m. For the 


case of n > 1 extra dimensions with small curvatures with the same value k, using results from 
the £> 0 (H 2 ) analysis one obtains constraints on the inverse of the curvature 1/k depending on n. 


For n = 2, 1/k < 1.5 X 10 +4 m; for n = 3, 1/k < 2.6 x 10 1 m; for n = 4, 1/k < 1.0 X 10 3 m; 


for n = 5, 1/Jfe < 3.9 X 10 

m. 
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m; for n = 6, 1/k < 4.3 x 10 6 m and for n = 7, 1/k < 8.8 X 10 


-7 


5. Comparison with other constraints 

The constraints obtained from molecular systems probe length scales in the order of Angstroms. 
This complements bounds probing subatomic to astronomical length scales obtained from other 








studies using distinct methodologies. The micrometer to millimeter range are probed in torsion- 
balance type experiments, with the tightest constraint obtained by Kapner et al. |20j for a single 
extra dimension of R < 4.4 x 10 5 m. 

Precision spectroscopies of hydrogen and muonic atoms have been interpreted in terms of 
the ADD-model [2lj resulting in typical constraints of R < 10 - ' 5 m. The interpretation is not 
straightforward because of the proton size puzzle [221 [23]; in fact, the argument has been turned 
around, and the existence of extra dimensions are instead invoked as a possible solution to the 
puzzle [2T]. In addition, unlike for molecules, in the treatment of atoms some assumptions had 
to be made on the wave function density at r = 0, typical for the s-states involved, causing 
problems in calculating the second integral of Eq. ([5]) over the electronic wave function that has 
a significant amplitude at r = 0 in atoms. 

To probe length scales in the subatomic range, one is ultimately limited by the increasing 
contributions from nuclear structure and the strong interaction. In contrast to QED calculations, 
the most accurate lattice-QCD calculation of light hadron masses only achieves relative 
accuracies in the order of a few percent. 

In high-energy particle collisions, higher-dimensional massive gravitons may be produced 
that could escape into the bulk, leading to events with missing energy in (3+l)-dimensional 
spacetime |3j, 25]. Based on this premise the phenomenology of the SN 1987A supernova was 
investigated, obtaining limits on extra dimensions of R < 3 X 10 6 m for n = 2, R < 4 X 10 7 
m for n = 3 and R < 2 x 10~ 8 m for n = 4 [ 26] . Similarly, from a missing energy analysis 
of proton colliding events at the LHC, the constraint R < 3.2 x 10~ 4 m can be extracted for 
n = 2 imposing the M4 +n = 1.93 TeV bound [27]. Also for n > 2, the bounds derived from the 
LHC are more stringent than those from molecules. However, additional assumptions beyond 
the ADD potential in Eq. ([3]) are necessary in order to interpret the LHC missing energy signals. 
Such assumptions are not needed for the molecular physics bounds, which are not sensitive to 
physics at very short distances. 

The higher-dimensional Planck mass M^ +n ^ and corresponding Planck length -Rpiyi+n) 
derived from the present constraints on R are also tabulated in Table [l] Note that in all cases the 
Planck energy scale in (4 + n) dimensions turns out to be in the range 1-500 MeV. Assuming a 
straightforward extrapolation of Eq. i> from molecular distance scales to nuclear distance scales, 
for an effect to be observable at molecular physics scales, the gravitational interaction must grow 
so rapidly at smaller distances that gravity becomes strong and non-perturbative already in the 
MeV scale. One would generally expect that strong gravity effects would already have been 
observed in scattering experiments, either through black hole creation, missing energy due to 
particles moving into the extra dimensions or other strong gravity effects, or large contributions 
to binding energies of nucleons and quarks. However, these expectations depend on a theory 
of gravity in the strong regime, and furthermore hadronic and nuclear binding is not a high- 
precision science, unlike atomic and molecular physics. Here we take the attitude that we are 
constraining the possibility that the gravitational force has a higher-dimensional dependence on 
r in the weak regime, where the computations are reliable, and we do not commit ourselves to 
what happens in the strong domain, where we cannot compute reliably. Our constraints are 
relevant only if gravity has a very soft behavior at shorter distances, so that its direct effects 
have escaped detection so far. 

The idea that gravity at short distances may be soft has been explored in four dimensions 
in several contexts. This includes asymptotically safe gravity [28], fat gravitons as a possible 
solution to the cosmological constant problem [29], a way of avoiding black holes and singularities 
in string theory [30] , and bouncing or self-inflating universes [31 . Therefore, given the limitations 
of our current understanding it is not unthinkable that this could happen in higher dimensions 
as well. 


6. Conclusion and outlook 

The alternative approaches for constraining compactification scales of extra dimensions, partially 
surveyed here, are all complementary as they probe different length and energy scales. Some 
approaches serve to produce tighter limits, however, often at the expense of additional 
assumptions. In the present study, the constraints are derived from precision measurements on 
molecules, leading to straightforward interpretations. Current state-of-the-art experiments on 
neutral hydrogen determine vibrational splittings of the order of ICC 4 cm' 1 , or 3 MHz [18] . Since 
the lifetimes of ro-vibrational quantum states in H 2 are of the order of 10 6 s [32], measurements 
of vibrational splittings of the order of 10 14 Hz could in principle be possible at more than 
20-digit precision, which leaves room for improvement “at the bottom” of over 10 orders of 
magnitude, if experimental techniques can be developed accordingly. 
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